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Abstract 

We propose a way of extending the Broad Histogram Monte Carlo 
method (BHMC) to systems with continuous degrees of freedom, and 
we apply these ideas to investigate the three-dimensional XY-model. 
Our method gives results in excellent agreement with Metropolis and 
Histogram Monte Carlo simulations and calculates for the whole tem- 
perature range 1.2 < T < 4.7 using only 2 times more computer effort 
than the Histogram method for the range 2.1 < T < 2.2. Our way 
of treatment is general, it can also be applied to other systems with 
continuous degrees of freedom. 

Usually one wants to calculate average values < Q >t of a certain quan- 
tity Q for a system in equilibrium at temperature T according to, for in- 
stance, the canonical ensemble 

_ Ee<Q>e 9(E) exp (-E/k B T) 
V E E 9(E)e W (-E/k B T) ' {±) 

where g(E) is the number of states with energy E, < Q >e denotes the 
micro-canonical average of Q at energy E, and k B is the Boltzmann constant 
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(in the rest ks = 1 is taken). It is clear that given g(E) and < Q >e one 
can calculate < Q >t at any desired temperature. 

The Broad Histogram Monte Carlo Method (BHMC) developed by de 
Oliveira et. al. Jl], ||] calculates g(E) and < Q >e directly. It chooses a 
microreversible protocol of allowed virtual movements in the space of states 
of the system, and counts N up (X) (Ndn(X)) as the number of allowed move- 
ments that increases (decreases) the energy of the configuration X by a fixed 
amount AE^ X . These movements are virtual, in the sense that they are never 
performed. They are introduced only to estimate g(E) and should not be 
mixed with the dynamics employed to take the samples. It has been shown 
H that the total number of possible ways to go up, summed over all states 
with energy E, equals the total number of ways to go down, summed over 
all states with energy E + AE^ X , i.e. 

g(E) < N up (E) >= g(E + AE Rx ) < N dn (E + AE &X ) > , (2) 

where < N up (E) > (< Nd n (E) >) is the micro-canonical average of N up (X) 
(N dn (X)) at energy E. By taking logarithms of both sides of Eq. (|2|) and by 
subsequent dividing by AE^ X one can obtain expressions for In g{E+AE^ x ) — 
In g(E) and (3{E) = din g{E)/dE that can be either added up or integrated 
numerically to obtain In g(E). The micro-canonical averages involved in the 
method can be calculated either analytically or by Monte Carlo simulations. 
In this case four histograms are required, i.e. N up (E), N^ n (E), Q(E) and 
the number of visits V(E). This procedure gives accurate results with small 
computer efforts for many discrete systems (Tj. ||. 

To extend the BHMC method to continuous systems we propose to 
choose a protocol of virtual random movements in the space of states of 
the system, such that for each allowed movement the probability to perform 
it equals the probability to revert it, i.e. 

-P(^old — ► ^ncw) = -P(^new _ ► ^old) ■ (3) 

Next, let us define the probability density function (p.d.f.) f x (AE) i.e. the 
probability to obtain an energy change between AE and AE + dAE starting 
from the configuration X. We propose 

N up (X) = f x (AE flx ) ; N dn (X) = f x (-AE fix ) . (4) 

Eq. (g) is still valid, and the method proceeds as before Q. 
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To test our ideas we chose the XY-model |J with the Hamiltonian 



n = - j 53 Si ■ Sj = - j 53 cos(e 4 - Oj) , (5) 

<y> <ij> 



where the summations < ij > are taken over all pairs of nearest-neighbor 
sites, 9i denotes the planar angle of the spin at site i and J is the maximal 
energy per bond (in the rest J = 1). As protocol of virtual random move- 
ments, we take at random one site i with angle 0i o id an d choose a new angle 
#inew for it with uniform probability on [— 7r,+7r]. We obtain Q 

/ x (a £ >4E/W) ; f(A*) = { Wi+A*) 2 |AS+ " I<A! ((5) 



otherwise 



where A denotes the number of sites i, and 



e< = - 53 ' 

i 



cos (0 iold - 9j) ; ^sW^«»ei) a +(E i sto9 i) 2 ( 7 ) 



The summation over j is always performed on all nearest neighbors of site 
i. Finally N up (X) (N dn (X)) equals the average of f(AE^) (f(-AB fc )) 
on all sites. 

To estimate the micro-canonical averages we binned the whole energy 
range of positive temperatures (i.e. energies per bond between —1.0 and 
0.0) in windows of size AEfi x and performed a micro-canonical simulation 
on each window as follows 0. Starting from a configuration inside the 
window one spin was rotated to a new angle chosen at random with uniform 
probability over the interval [— n, ir]. If the new system energy falls inside 
the window, the change is performed, otherwise it is rejected. It is clear 
that this sampling process maintains a detailed balance condition with equal 
probabilities for the old an the new configurations and, therefore, it samples 
all the states inside the window with the same probability. An entire lattice 
sweep was performed by repeating this procedure on all sites. 

Fig. Q compares the specific heat we obtained for the SC 10 x 10 x 10 
XY-model with periodic boundary conditions by using the BHMC method 
(circles), the Histogram Monte Carlo method (HMC) |?], ^] (squares) and 
Metropolis || simulations (diamonds). For the BHMC method A.Eu x = 6.0 
(AEfi x = 0.002 in energy-per-bond units) was chosen. 

The HMC method took 188 sec. per run, a Metropolis for one temper- 
ature point 28.5 sec. and the BHMC method 383 sec. per run on a Digital 
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Figure 1: Specific heat for the SC 10 x 10 x 10 XY- model obtained from the BHMC 
method (circles), the HMC method (small squares) and Metropolis simulations (filled dia- 
monds). The error bars on the left (right) side correspond to 1.0 (3.5, i.e. 99% confidence 
level) standard deviations for eight runs. 



Alpha workstation, i.e. only 2.0 times more computer effort than the HMC 
method and less than 14 Metropolis simulations. It can be observed that all 
methods coincide on the temperature range 2.1 < T < 2.2, i.e. the classical 
range of validity for this HMC simulation [1C], but the BHMC method is 



still precise on the whole temperature range 1.2 < T < 4.7. 

In conclusion, our proposal to apply the BHMC method to continuous 
systems gives excellent agreement with Metropolis and HMC simulations. 
The strategy proposed is completely general, it can also be applied to other 
systems with continuous degrees of freedom. 



References 

[1] P.M.C. de Oliveira, T.J. P. Penna and H.J. Herrmann, Braz. J. of Phys. 
26, 677-683 (1996); also in |cond-mat/9610041| . 

[2] P.M.C. de Oliveira, Int. I. Mod. Phys. C9, 497 (1998). 



[3] P.M.C. de Oliveira, to appear in Eur. Phys. J. B; also in cond- 
mat/ 9807354. 



[4] J.D. Muhoz and H.J. Herrmann, to appear in Int. J. Mod. Phys. C; 
also in |cond-mat/9810024 . 



[5] J. Adler et. al, Physica A 201, 581 (1993). 



4 



[6] N. Metropolis et. al., J. Chem. Phys. 21, 1087 (1953). 

[7] Z.W. Salsburg et. al, J. Chem. Phys. 30, 65 (1959). 

[8] A.M. Ferrenberg and R.H. Swendsen, Phys. Rev. Lett. 61, 2635 (1988). 

[9] W.H. Press et. al., Numerical Recipes in Fortran 90, 2nd. Ed. (Cam- 
bridge University Press, New York, 1996). 

[10] M.E.J. Newman and R.G. Palmer, Error estimation in the histogram 
Monte Carlo method, Preprint cond-mat / 9804306 . 



5 



